AN  ADAPTIVE 

KALMAN  FILTER  EXCISOR  FOR 
SUPPRESSING  NARROWBAND  INTERFERENCE 


,  dtic 

CLECTE 
'•X  i- £3  2  8 1994 

?>  ri 


by 

Brian  W.  Kozminchuk 


This  doj’22i=n»  hiS  re«a  Qpptov^d 
lot  publir  isisi5e  cad  jol*.  its 
diitnbutioa  is  iiauaiiied. 


9±Q6299 

ilii 


DEFENCE  RESEARCH  ESTABLISHMENT  OTTAWA 

TECHNICAL  NOTE  93-20 


Canada 


rrc  C' 


a:  .-D  i 


December  1993 
Ottawa 


n 


o=; 


00 


3  U 


Best 

Available 

Copy 


m  National 
Oalanc* 


Daiansa 

nationjia 


AN  ADAPTIVE 

KALMAN  FILTER  EXCISOR  FOR 
SUPPRESSING  NARROWBAND  INTERFERENCE 


by 


Brian  W.  Kozminchuk 

Electronic  Support  Measures  Section 
Electronic  Warfare  Division 


DEFENCE  RESEARCH  ESTABLISHMENT  OTTAWA 

TECHNICAL  NOTE  93-20 


PCN 

041LK11 


December  1993 
Ottawa 


ABSTRACT 


Severed  reports  have  been  written  by  the  author  which  characterize  the  performance  of 
a  Kedman  filtering  technique  to  suppress  narrowband  interference  from  direct  sequence 
spread  spectrum  communications  systems  in  a  non-adaptive  context,  i.e.,  when  the  user 
has  some  a  priori  knowledge  of  the  interference  characteristics.  This  report  expands  on 
this  technique  by  presenting  an  adaptive  scheme  which  is  useful  for  the  situation  when 
the  interference  characteristics  are  unknown.  In  this  context,  the  Kalman  filter  must 
“learn”  to  achieve  optimal  performcince  through  the  adjustment  of  one  of  its  parame¬ 
ters.  The  criterion  for  optimality  is  the  minimization  of  the  mean-squared  error  at  the 
output  of  the  canceller,  where  this  error  consists  of  spread  spectrum  signal,  noise  and 
residual  interference.  The  reasonable  assumption  is  made  that  minimizing  the  mean  of 
this  squared  error  with  respect  to  the  appropriate  Kalman  filter  parameter  is  equivalent  to 
minimizing  the  mean  squared  value  of  the  residual  interference.  Examples  of  the  dynamic 
behaviour  of  the  adaptive  interference  suppressor  are  presented  for  narrowband  Gaussian 
noise  interference  with  bandwidths  ranging  from  1%  to  5%  of  the  chip  rate. 

RESUME 

L'auteur  de  ce  rapport  a  deja  public  plusieurs  documents  sur  I’utilisation  de  filtres  de 
Kalman  pour  annuler  les  I’interferences  a  bandes  passantes  etroites  dans  des  systemes 
de  communication  utilisant  des  signaux  a  spectres  etales  par  sequences  directes.  Le 
cadre  de  ces  travaux  etait  liinite  a  des  techniques  non  adaptatives  oil  les  caracteristiques 
de  I’interference  etaient  connues,  a  priori,  de  I’utilisateur.  Ce  rapport  presente  des 
developpements  oil  une  approche  adaptive  est  utilisee  pour  annuler  des  interferences  in- 
connues.  Le  filtre  de  Kalman  doit  alors  “apprendre”  a  ajuster  un  de  ses  parametres  pour 
effectuer  le  meilleur  trutement.  L’erreur  est  constituw  du  signzJ  a  spectre  etale,  du  bruit 
et  de  I’interference  residuelle,  et  la  minimisation  de  I’erreur  quadratique  moyenne  a  la 
sortie  de  Tannuleur  sert  de  critere  d’optimisation.  On  fait  rhypoth«e  que  de  minimiser 
I’erreur  quadratique  moyenne  par  rapport  au  parametre  approprie  du  filtre  de  Kalman 
^uivaut  a  minimiser  la  valeur  quadratique  moyenne  de  I’interference  residuelle.  Des 
exemples  de  la  dynamique  du  comportement  de  I’annuleur  adaptatif  d’interference  sont 
presentes  pour  de  I’interference  gaussienne  a  bande  passante  etroite  (1%  a  5%  du  debit 
numerique). 
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EXECUTIVE  SUMMARY 


This  technical  note  presents  an  adaptive  Kalman  filtering  technique  to  suppre-ss  iinirow- 
band  interference  from  direct  sequence  spread  spectrum  communications  systems.  This 
is  an  expansion  of  the  non-adaptive  Kalman  filtering  approach  reported  on  by  the  author 
in  earlier  reports. 

Spread  spectrum  signals  are  used  extensively  in  military  communication  sys¬ 
tems.  The  technique  described  herein  applies  equally  to  both  Electronic  Support  .Mea¬ 
sures  (ESM)  systems  and  direct  sequence  spread  spectrum  communication  systems.  In  the 
former  application,  the  ESM  system  may  be  attempting  to  intercept  the  spread  spectrum 
signal,  but  the  narrowband  interference  may  be  hampering  this  effort.  In  the  latter  appli¬ 
cation,  the  spread  spectrum  communication  system  may  require  additional  assistance  to 
suppress  the  interference.  Since  the  open  literature  has  been  devoted  to  this  latter  case, 
the  material  presented  here  focuses  on  this  application. 

One  of  the  attributes  of  direct  sequence  spread  spectrum  communication  systems 
is  their  ability  to  combat  interference  or  intentional  jamming  by  virtue  of  the  system's 
processing  gain  inherent  in  the  spreading  and  despreading  process.  The  interference  can 
be  attenuated  by  a  factor  up  to  this  processing  gain.  In  some  cases  the  gain  is  insufficient 
to  effectively  suppress  the  interferer,  leading  to  a  significant  degradation  in  system  per¬ 
formance  as  manifested  by  a  sudden  increase  in  bit  error  rate.  If  the  ratio  of  interference 
bandwidth  to  spread  spectrum  bandwidth  is  small,  the  interference  can  be  filtered  out  to 
enhance  system  performance.  However,  this  is  at  the  expense  of  introducing  some  dis¬ 
tortion  onto  the  signal.  This  process  of  filtering  is  sometimes  referred  to  as  interference 
excision. 

The  adaptive  technique  reported  upon  herein  is  useful  for  the  situation  when  the 
interference  characteristics  are  unknown.  In  this  context,  the  Kalman  filter  must  “learn” 
to  achieve  optimal  performamce  through  the  adjustment  of  one  of  its  parameters.  The 
criterion  for  optimality  is  the  minimization  of  the  mean-squared  error  at  the  output  of  the 
canceller,  where  this  error  consists  of  spread  spectrum  signal,  noise  and  residual  interfer¬ 
ence.  The  reasonable  assumption  is  made  that  minimizing  the  mean  of  this  squared  error 
with  respect  to  the  appropriate  Kalman  filter  parameter  is  equivadent  to  minimizing  the 
mean  squared  value  of  the  residual  interference.  Examples  of  the  dynamic  behaviour  of  the 
adaptive  interference  suppressor  are  presented  for  narrowband  Gaussian  noise  interference 
with  bandwidths  ranging  from  1  to  5%  of  the  chip  rate. 
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1.0  INTRODUCTION 


Direct  sequence  spread  spectrum  rommunic.ation  systems  have  an  inher<Mit  proressing 
gain  which  can  reduce  the  effects  of  jammers  or  unintentional  interference.  When  t  In  sc 
intruding  signals  have  a  power  advantage  over  the  spread  spectrum  system,  a  severe  degra¬ 
dation  in  communications  results.  However,  communications  can  be  enhanced  somewhat 
by  filtering  the  interference,  particularly  if  its  bandwidth  is  significantly  less  than  the 
bandwidth  of  the  spread  spectrum  signal. 

Several  reports  have  been  written  by  the  author  [1,  2,  .3,  4,  .5]  on  the  performance 
characteristics  of  a  Kalman  filtering  approach  to  narrowband  interference  suppression  in  a 
non-adaptlve  context,  i.e.,  when  the  user  has  some  a  priori  knowledge  of  the  interference 
characteristics.  This  report  focusses  on  an  adaptive  scheme  which  is  useful  for  the  situation 
when  these  characteristics  are  unknown,  i.e.,  through  some  criterion,  the  Kalman  filter 
must  adjust  its  parameters  to  achieve  optimal  performance. 

The  outline  of  the  report  is  as  follows.  Sections  2.0  and  3.0  review  the  basic 
elements  of  the  non-adaptive  interference  suppression  system,  with  Section  2.0  describing 
the  spread  spectrum  communications  model  and  Section  3.0  detailing  the  interference 
estimator.  This  latter  section  includes:  a  description  of  the  state  space  model  which 
generates  a  function  related  to  the  interference;  the  extended  Kalman  filter  algorithm 
and  resultant  digital  phase-locked  loop  (DPLL)  which  is  used  to  estimate  the  interferer’s 
phase;  and  the  envelope  estimator  which  is  combined  with  the  output  of  the  Kalman  filter 
to  produce  an  estimate  of  the  interference.  Section  4.0  describes  the  adaptive  algorithm 
which  is  based  on  minimizing  the  mean-squared  value  of  the  error  signal  at  the  output  of 
the  canceller  while  one  of  the  Kalman  filter  parameters  is  adjusted  over  time  in  a  controlled 
fashion.  This  secticr.  also  includes  ciinulation  results  for  the  case  of  narrowband  Gaussian 
noise  of  different  bandwidths  corrupting  the  spread  spectrum  signal.  Finally,  Section  5.0 
summarizes  the  results  of  this  report  and  suggests  areas  for  further  research. 

2.0  COMMUNICATIONS  MODEL 

The  basic  elements  of  the  BPSK  PN  spread  spectrum  system  are  shown  in  Fig.  1.  The 
re  eived  waveform  r{t),  consisting  of  a  spread  spectrum  signal,  additive  white  Gaussian 
noise,  and  narrowband  interference  is  applied  to  a  bandpass  filter  with  the  transfer  func¬ 
tion  Hhp(f),  whose  output  is  defined  as 


u(t)  s  s(t) -b  n(f) -I- 1(0. 


(1) 


PN 

Cortelator 


Figure  1:  Spread  spectrum  communications  model. 

The  bandpass  filter  for  the  application  considered  here,  is  assumed  to  be  a  filter 

matched  to  a  chip  and  centered  at  the  carrier  frequency  jJq  of  the  spread  spectrum  signal 
defined  as 

3(t)  =  a{t)  cos(u^t)  (2) 

where 

=  (3) 

* 

In  Eq.  (3),  Di,  is  a  sequence  of  data  bits  of  amplitude  (±1),  Ti,  is  the  reciprocal  of  the  bit 
rate  R^,  and  bk{t  —  kT^)  is  the  PN  sequence  pattern  for  the  bit,  i.e., 

bk{t)  =  Zc,Mt-jTc)  (4) 

with  L  being  the  number  of  pseudo  random  chips  per  bit,  or  the  processing  gain,  ctj  is 
the  code  sequence  for  the  bit,  Te  is  the  recipro  :al  of  the  chip  rate  Rc  and  q{t)  represents 
the  basic  chip  pulse  of  energy  £’«. 

The  noise  n(t)  in  Eq.  (1)  is  Gaussian  and  has  a  power  spectral  density 

5n(/)  =  Y|/f6p(/)l*,  10) 

where  Nq/2  is  the  power  spectral  density  of  the  assumed  white  Gaussian  noise  from  the 
channel.  The  interference  term  in  Eq.  (1)  is  of  bandwidth  Bj  <<  2Rc  and  is  defined  as 

j(t)  =  I{t)  cos{uot  -hff(t))  (6) 


*> 


vlnTt’  l{()  i«  tiu-  mlfrfcnTicf' aikI  0{t)  is  (In-  ph.isc  iruxlul.ii ion.  It  )i,is  Ix-rn 
•ussninoil  t  iiAt  I  h<'  t‘ir<'ci  of  tin*  liAntIp.uss  (ilt«‘r  //•,p(  J  )  i*-  no”  I  ini  Me  on  i  li<*  inti-rf<‘r<-n<'-  lit). 

Roforrini'  to  Fit'.  1.  llio  output  ii(/)  of  tli<*  o.unlii.i's  tiltor  //»,,,(/)  of  nAin  y  l.^T- 
IS  bandpass  sampU'd  and  applied  to  a  limiter/hamipass  filter  and  intcrfercnic  estimator. 

Consider  the  bandpass  sampler  first.  The  analog  signal  n(t)  from  Eq.  (1)  is 
sampled  at  /,  =  2Rc  (m/,  =  wo/‘2jr  +  for  some  integer  m).  The  resultant  sampled 
signal  is,  therefore, 

tin  =  -Jn  +  +  »n.  (f) 

where  s,  consists  of  the  sequence  {. . .  ,0,  (  — 1  )"nn,0.  (  -  1  . . .)  where  the  a,  are 

of  energy  Eg  and  co<led  according  to  CkjDk  for  the  j'*'  chip  in  the  transmitted  bit, 
n„  *  arc  uncorrelated  Gaussian  noise  samples  of  variance  =  £‘c(.Vo/2),  and  in  is  the 
sampled  version  of  Eq.  (6).  The  samples  u,  are  applied  to  the  interference  estimator  and 
interference  canceller. 

Consider  now  the  branch  containing  the  limiter.  Here,  Un  is  applied  to  a  lim¬ 
iter/bandpass  filter.  The  input  to  the  limiter  referenced  to  the  interference  is  redefined 

as  _ 

Un  =  \/(/n  -f  n'l.n  +  Ul.nl’  +  Kj.n  +  +  «,  +  Ou.n)  (8) 

where 

<Pu,„  =  arctan 

is  a  noise-like  phase  fluctuation  on  the  interferer’s  phase  and  is  due  to  the  noise  and 
spread  spectrum  signal.  The  terms  n',  a',  „  and  „  are  in-phase  and  quadrature 

components  of  the  noise  and  spread  spectrum  signal  with  respect  to  the  interference  phase 
6n-  The  output  of  the  limitrr/bandpass  filter  is  (6) 

4/4' 

6pn  =  —  cos(w>on  +  6„  +  c„,„],  ( 10) 

where  A'  is  the  limiter’s  output  level.  This  signal  is  redefined  as  (letting  A'  =  >/2t/4  for 
convenience) 

bpn  =  ^'2  cos[a;on  -f  4-  0„,„) .  (11) 

It  should  be  noted  that  for  large  interference-to-noise  ratios  in  which  the  interference  is 
of  constant  envelope,  in  Eq.  (11)  is  approximately  Gaussian  [6].  The  sampled  signal 

'Coherent  bandpass  sampling  has  been  assumed,  so  that  the  in-phase  component  is 
(ni.n/\/2)cos(nr/2)  and  the  quadrature  component  is  (nj,„/\/2)sin(nT/2> 


-t-  .1',  „ 

+  Ul.n/ '/5  +  rt'i  „ 


(9) 


.3 


in  Eq.  (ll)  is  what  is  processed  by  the  Kalman  fitter,  which  estimates  the  phase  0^  of  the 
interference. 

3.0  THE  INTERFERENCE  ESTIMATOR 

This  section  presents  the  interference  estimator.  It  consists  of  the  extended  Kalman 
filter  which  estimates  the  phase  of  the  interference,  and  the  envelope  estimator  which  is 
combined  with  the  Kalman  filter  result  to  produce  an  estimate  of  the  interference. 

First,  the  analog  and  digital  state  space  models  on  which  the  Kalman  filter  algo¬ 
rithm  is  based  will  be  presented.  This  will  be  followed  by  the  extended  Kalman  algorithm 
used  to  estimate  the  state  variables.  Finally,  the  section  concludes  with  two  structures  for 
estimating  the  interference  using  the  results  of  the  Kalman  filter  and  envelope  estimator. 


3.1  THE  STATE  SPACE  MODEL 


The  development  here  follows  that  in  (7).  A  block  diagram  of  the  analog  state  space  model 
which  generates  the  analog  version  of  a  function  related  to  in  Eq.  (11)  is  shown  in 
Fig.  2;  this  function  is  represented  by  bp{ty  ^  which  is  defined  later  in  this  section.  The 
state  space  equations  are  defined  as 


ii(0 

ij(<) 


0  d 

0  —a 


xi(<)  , 


w(t). 


(12) 


* 


where  xi{t)  and  X2(f)  are  the  phase  function  and  modulating  signed  state  variables  of  the 
interference,  respectively.  In  compact  form  Eq.  (12)  becomes 

i(f)  =  Ax(t)  +  (13) 

The  interferer’s  modulating  signal  Z3(()  is  assumed  to  be  generated  by  applying  Gaussian 
white  noise,  w{t),  of  unit  variance  to  .an  amplifier  with  gain  y/KZ  followed  by  a  first  order 
low  pass  filter  of  bandwidth  a  =  2tq/.  The  steady  state  variance  of  the  coloured  noise 
at  the  output  of  the  lowpass  filter  is  unity  if  —  2a.  The  output  is  multiplied  by  the 
frequency  deviation  constant,  d,  and  then  integrated  to  yield  the  phase  function 

xi(t)  =  df  X2(T)dr.  (14) 

Jo 

*In  the  analog  version  of  Eq.  (11),  the  noise  is  one  of  the  arguments  of  the  cosine  function.  For  the 
model  discussed  here,  the  noise  will  be  additive. 
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An  FM  type  of  signal  (interfcrer)  is  generated  by  plias(>  modiilaling  a  carrier,  -.o,  which 
has  been  set  equal  to  the  carrier  of  the  spread  spectrum  signal,  since  any  olfsets  from  -o 
can  be  incorporated  into  Xi(<).  yielding 

j/(f)  =  V2cos(u;ot  +  xi(0l-  (b’j) 

This  signal  is  corrupted  by  white  Gaussian  noise  v{t)  with  power  spectral  density  N„i„/2 
generally  different  from  Nq/2  defined  earlier  for  the  communications  model  in  Fig.  1.  The 
result  is  an  FM  signal  related  to  the  interference  in  Eq.  (11)  except  for  the  fact  that  the 
noise  term  is  additive  here,  i.e.. 


bp(t)  =  v(t).  (16) 

For  large  interference- to-signal-plus-noise  conditions,  the  final  effect  of  the  additive  noise 
in  ^timating  the  phase  xi(0  in  Eq.  (15)  will  be  similar  to  the  effect  of  <i>u.n  in  Eq.  (11). 

Given  the  observation  process  bp*(t)  from  Eq.  (16),  the  objective  is  to  estimate 
the  state  ii(t)  and  the  related  state  ia(0,  as  well  as  y[t)  in  Eq.  (15).  One  approximate 
technique  reserved  for  nonlinear  estimation  problems  [S]  is  the  extended  Kalman  filter 
algorithm. 


v(0 


Figure  2:  Block  diagram  of  the  state  space  model. 

Before  the  algorithm  is  presented,  the  discrete  forms  of  Eqs.  (12)  and  (13)  will 
be  stated,  namely, 

Xn+t  =  +  w„,  (17) 

where  n  represents  discrete  time, 

i  Oil 
0  Oaa 


1  ^(l-e-"^*) 

0  e-“^* 
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(18) 


u>n  =  /  <j>gw(z)dz. 

JnT. 


In  Eq.  (18),  0  =  d/a  and  is  termed  the  bandwidth  expansion  ratio  in  units  of  volts'^ 
[9j.  In  Eq.  (19),  ti)n  is  a  stationary  zero-mean  white  Gaussian  vector  sequence  whose 
covariamce,  (i.e.,  £{iiJBto|,}),  is  given  by 


r{n+t)T, 

I  <j>gg*<j>du. 
JnT, 


Finally,  Eq.  (16)  becomes. 


¥n  =  y«  +  Vn> 


where  y„  is  the  sampled  version  of  Eq.  (15),  i.e., 

j/n  =  >/2  cos(a?on  +  ii,n)  (22) 

It  has  been  assumed  that  bp'{t)  from  Eq.  (16)  has  been  applied  to  a  bandpass  filter  of 
bandwidth  B,  whose  output  has  been  sampled  at  a  rate  so  as  to  yield  uncorrelated  noise 
samples  Vn  of  variance  al  =  B(Noi,/2)  (10). 

The  noise  covariance  matrix,  is  obtained  by  substituting  Eq.  (18)  and  g 
from  Eq.  (13)  into  Eq.  (20)  and  integrating  each  term,  yielding 

V  /?  (l 1 


^  (l  - 


1  _  e~**h 


where  7  =  2v/{aT,),  which  is  the  ratio  of  the  sampling  rate  /,  =  1/T,  to  the  interferer’s 
modulating  signal  bandwidth  a/.  The  discrete  form  of  the  interference  model  is  shown  in 
Fig.  3. 


3.2  EXTENDED  KALMAN  FILTER  EQUATIONS 

The  extended  Kalman  filter  equations  [8j  can  now  be  applied  to  the  discrete  model.  The 
resulting  algorithm  is  listed  in  Table  1.  In  the  algorithm  in  is  defined  as  the  estimate  of 
the  true  state  vector  x„,  in  is  defined  as  the  post  error  x„  —  x„,  Kn  is  the  Kalman  gain 
vector,  Vf,  =  £{x„x^}  is  the  post  error  covariance  and  =  E{x„|„_ix|,|„_,}  is 

the  prior  error  covariance.  Furthermore,  fZ\n-\  **  known  as  the  “coupled  equation",  since 
the  estimate  of  the  interferer’s  phase  is  contained  in  it.  However,  this  equation  can  be 
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Figure  3:  Discrete  form  of  the  interference  model. 


further  approximated  as  noted  in  [7].  The  approximation  is 


3.3  THE  DIGITAL  PHASE-LOCKED  LOOP  (DPLL) 

The  Kalman  algorithm  results  in  a  digital  PLL  {7|.  This  fact  will  be  elaborated  upon 
here  so  as  to  illustrate  the  sampling  rate  requirements  which  are  related  to  the  allowable 
interference  bandwidths  which  the  system  is  capable  of  handling. 

Substitution  of  the  Kalman  gain  Kn  from  Table  (1)  into  the  state  estimate 
equation  for  in  yields 

in  =  «^*n-i  +  V  |-6pJ,N/2sin[u;on  + 

+  sin[wbn  +  ^n|n-ll  cos(won  +  ^n|n-l]}  (25) 

Consider  the  term  within  braces.  Substituting  for  6pJ,  from  Eq.  (21)  and  noting  that 
xjn  =  one  obtains  the  intermediate  variable 

kn  =  sin(0„  -  5„|n_i)  -  sin(2u;on  +  +  9n\n~i ) 

+  sin(2a;ori  +  20„)  -  >/2un  sin(u;on  +  )  (26) 

which  consists  of  four  terms.  The  first  refers  to  the  baseband  term.  The  second  and  third 
are  second  harmonic  terms.  Recalling  that  t7„  is  the  sampled  version  of  bandpaiss  noise 
centered  at  wo,  the  fourth  term  turns  out  to  be  baseband  noise  which  can  be  modeled  as 
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a  zero  mean  white  Gaussian  noise  process  (9j.  The  first  and  fourth  terms  are  the  desired 
ones,  since  they  correspond  to  their  analog  PLL  counterparts.  Also,  as  in  the  analog 
implementation,  the  second  harmonic  terms  can  be  neglected,  provided  they  do  not  alias 
into  the  baseband  region. 

To  ensure  that  no  spectral  aliasing  of  these  terms  occurs,  the  sampling  rate  must 
be  selected  according  to  /,  >  4Bi.  This  can  be  seen  from  the  following  heuristic  argument. 
As  a  starting  point  consider  analog  spectra.  Then  assume  that  no  noise  is  present  and 
^(t)  St  9{t).  Also  let  the  terms  within  braces  in  Eq.  (25)  be  represented  as  follows: 


9i(<)  =  6p'(f)  =  v'2cos[wot  +  5(01 

(27) 

92(f)  =  N/2sin[wot  +  5(t)] 

(2«) 

93(f)  =s  sin[wof  +  5(1)]  cos[wot  -i-  5(t)) 

=  -sin[2a;of +  25(t)). 

(29) 

The  amplitude  spectrum  of  qi{t)  is  shown  in  Fig.  4(a).  It  is  assumed  that  the  amplitude 
spectrum  of  qj(t)  is  approximately  equal  to  it  for  the  case  of  9{t)  ss  0(t).  Taking  the 
product  of  gi(t)  and  92(f)  in  (2T)  and  (28)  results  in 

94(1)  =  sin(d(t)  —  0(t))  +  sin[2wof  +  ^(f )  +  )]  (30) 


whose  spectrum  is  shown  in  Fig.  4(b),  the  convolution  of  Qi(/)  2uid  Q2{f)-  If  an  integer 
multiple  of  the  sampling  frequency  mf,  is  chosen  as  shown  in  Fig.  4(b)  then  the  required 
sampling  rate  for  no  aliasing  in  the  baseband  region  is  /,  >  4Bi.  Thus,  the  baseband  and 
second  harmonic  terms  interleave  as  shown  in  Fig.  4(c).  Similarly,  the  spectrum  for  93(4) 
and  its  sampled  spectrum  are  shown  in  Figs.  5(a)  and  (b). 

The  second  harmonic  terms  can  therefore  be  neglected  and  Eq.  (25)  becomes 


0x„_}  +  V 


0 


+  Vi,n 


{sin[^  -  +  «;,]} 


(31) 


where 

Vn  =  -v/2w„sin(w(,n  + 


(32) 
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Flrequencf  Qiz> 

(•) 


(b) 


Figure  4:  lUustratioa  of  the  sampling  rate  requirements:  (a)  Hypothetical  interference 
spectrum  of  9i(0i  (b)  Convolution  of  Q\{f)  with  Q2{f);  (c)  Sampled  spectrum  of  the 
product  q\{t)q2{t). 
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Figure  5:  Analog  sind  sampled  spectra  of  93(f)  (see  Eq.  (29)). 
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and  is  of  variance  {N^,f2)f,. 

The  state  space  diagram  of  the  DPLL  is  illustrated  in  Fig.  6.  Observe  that 


Figure  6:  State  space  representation  of  the  DPLL. 

the  phase  estimate  is  applied  to  a  phase  modulator  of  frequency  wo,  resulting  in  a 
“VCO”  signal  defined  as 


VCOn  =  ->/28in(won  +  ^n|n-l)- 


(33) 


Shifting  this  signal  by  —1(12  radians  results  in  a  filtered  estimate  of  in  Eq.  (21). 

It  should  be  noted  that  the  Kalman  algorithm  in  Table  (1)  which  is  based  on 
the  state  space  model  shown  in  Fig.  2  will  le£ui  to  a  steady  state  error  if  the  phase  is  a 
ramp.  The  steady  state  error  can  be  reduced  to  zero  by  modifying  the  interference  model 
[11].  The  modification  involves  replacing  the  low  pass  filter  in  Fig.  2  by  an  integrator 
(i.e.,  by  setting  a  =  0  in  A  in  Eq.  (13)).  With  this  change,  the  state  transition  matrix  <i> 
and  covariance  matrix  V„  become,  respectively. 


0  = 


'  1  dr/ 

_  1  <l>12 

0  1 

L  0  (^2 

(34) 


(35) 


and 

y  20a^Tj 

2a  2/3a*r;  2aT, 

Using  these  terms  in  the  Kalman  equations  results  in  a  slightly  modified  algorithm.  This 
is  the  2dgorithm  that  has  been  used  in  the  simulations. 
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3.4  INTERFERENCE  ESTIMATOR  STRUCTURES 

This  section  presents  two  structures  for  the  interference  estimator  shown  in  Fig.  1  and 
which  includes  the  Kalman  filter  discussed  in  the  previous  section.  These  structures  have 
been  used  in  connection  with  the  non-adaptive  Kadman  algorithm  [5, 12]  and  are  presented 
here  for  later  reference. 

As  noted  earlier,  the  Kalman  algorithm  produces  a  filtered  estimate  of  t/n  in 
Eq.  (22),  which  can  be  used  as  a  basis  for  estimating  a  sampled  version  of  the  envelope 
of  t(t)  (i.e.,  In  in  Eq.  (6))  which  is  mixed  with  the  same  yn  to  yield  an  estimate  of 
the  interference.  The  first  structure  is  a  .straightforward  implementation  of  this  idea. 
The  second  structure  carries  out  additional  phase-smoothing  on  the  Kalman  filter  phase 
estimate;  this  adds  extra  complexity,  but  improves  the  perform2uice  somewhat  [12]. 

The  first  structure  is  illustrated  in  Fig.  7.  The  output  of  the  first  multiplier  is  a 
baseband  term  and  is,  using  Eq.  (7),  excluding  the  y/2  factor, 

/'n  =  /nCOs(5n-5n|i.-l) 

+  n„  cos(won  +  5„jn-i )  +  a*  cos(fl„|„>, ) 

+  anCos(2u^n-|-d„|„_i).  (36) 

Equation  (36)  consists  of  four  terms:  the  first  term  is  related  to  the  desired  amplitude  of 


Figure  7:  Block  diagram  of  the  interference  estimator. 

the  interference;  the  second  is  approximately  Gaussian  baseband  noise  [9];  and  the  third 
and  fourth  terms  are  noise-like  terms  emanating  from  the  spread  spectrum  signal.  The 
fourth  term,  because  of  the  sampling  rate  conditions  discussed  in  Section  3.3,  is  essentially 
filtered  out  by  the  low  pass  filter  of  bandwidth  Bcpf  <  0.50  Hz  and,  therefore,  will  be 
Ignored  in  the  baseband  simulations  to  be  discussed  in  the  next  chapter.  The  term  \/2 
is  filtered,  resulting  in  the  estimate  of  the  interference  envelope,  /„/\/2.  Combining  this 


i 
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with  an  estin  'te  of  y„  in  Eq.  (22),  y„,  shown  in  Fig.  7  yields  the  estimate  of  the  interference 

A  *  * 

in  =  /nCOsfwon  +  <?n|n-l),  (37) 

which  is  subtracted  from  Un  as  illustrated  in  Fig.  1.  Because  the  input  was  coherently 
bandpass  sampled  at  twice  the  chip  rate,  the  output  of  the  canceller  (summer)  undergoes 
decimation  and  alternate  sign  reversal,  resulting  in  the  error  sign2J  which  undergoes 
correlation  with  the  local  PN  code  and  bit  detection. 

The  second  structure  is  illustrated  in  Fig.  8.  There  were  two  reasons  for  intro¬ 
ducing  this  structure.  The  hrst  was  to  counter  the  effects  of  the  group  delay  produced  by 
the  Butterworth  low  pass  filter  in  Fig.  7  by  using  a  linear  phase  FIR  filter  with  known 
delay  D.  The  second  was  to  do  some  additional  smoothing  on  the  phase  estimate  dnin-i 
from  the  Kalman  Alter  which  would  introduce  the  same  delay  D.  With  phase-smoothing 
and  the  FIR  low  pass  Alter,  the  interference  estimate  becomes  tn-o  which  is  subtracted 
from  a  delayed  version  of  the  input,  Un-p.  The  result  is  the  error  signal  Cn-p  which  is 
applied  t't  the  PN  correlator  whose  PN  sequence  is  delayed  by  D  samples  as  well. 


Figure  8:  Interference  estimator  when  phase-smoothing  is  being  used. 


4.0  AN  ADAPTIVE  KALMAN  FILTER  EXCISOR 


To  this  point,  a  non-adaptive  Kalman  filter  has  been  presented.  It  will  be  shown  that 
by  varying  the  frequency  deviation  constant  d,  while  keeping  the  variance  of  the  ob¬ 
servation  noise  and  the  modulating  signal  bandwidth  a/  of  the  Kalman  filter  model  for 
FM  interference  fixed,  a  minimum  for  the  mean-squared  value  of  the  residual  interfer¬ 
ence,  can  be  obtained.  This  suggests  an  adaptive  scheme  which  this  section  is 

concerned  with. 

The  proposed  scheme  is  to  vary  d  over  time  (with  cr^  and  a/  fixed)  and  to  monitor 
an  estimate  of  EfeJ},  the  mean-squared  value  of  the  excisor’s  error  signal  in  Fig.  1,  which 
is  related  to  as  discussed  below. 


0  200  400  600  800  1000  1200  1400  1600  1800  2000 

Time  (sec.) 

Figure  9:  Profile  of  the  frequency  deviation  constant  d  used  in  the  phase-  and  envelope¬ 
tracking  tests. 


As  an  example,  Fig.  9  illustrates  a  time-varying  frequency  deviation  constant 
d.  The  interference  for  this  case  was  narrowband  Gaussian  noise  with  bandwidth  Bi  = 
0.01  Hz  and  the  interference  estimator  was  the  system  in  Fig.  7  with  Blpf  =  0.40  Hz. 
The  phase-  and  envelope-tracking  capabilities  of  a  combination  of  the  Kalman  filter  and 
envelope  estimator  arising  from  this  time- varying  function  are  shown  in  Figs.  10  and  11. 


Referring  to  Fig.  10,  for  large  values  of  d,  the  phase-noise  is  evident  until  at  about  400 
iterations,  at  which  point  the  phase-noise  has  decreased  substantially,  indicating  that  tlie 
bandwidth  of  the  Kalman  filter  has  decreased.  Observe  also  that  initially,  the  algorithm 
exhibits  some  slippage,  taking  from  10  to  15  radians  before  it  settles  down  to  a  tracking 
error  of  2t  radi2uis  after  approximately  50  iterations;  the  slippage  is,  in  part,  due  to  d 
being  initially  too  large.  At  the  other  extreme,  when  d  becomes  too  small,  the  algorithm 
eventually  loses  track  because  the  rate  of  change  of  the  phase  is  too  fast  for  it  (i.e.,  the 
bandwidth  of  the  Kalman  filter  is  too  small). 


Figure  10:  An  example  of  the  phase-tracking  capability  of  the  Kalman  filter  when  the 
frequency  deviation  constant  d  was  varied  according  to  the  profile  in  Fig.  9. 
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Figure  11  illustrates  the  estimate  of  the  envelope  of  the  interferem  e  <liirin;'  tli<- 
same  time  period.  As  d  decrea-ses,  !„  appears  to  improve,  until  at  approxiin.aiely  1200 
iterations,  the  estimate  degrades  substantially  due  to  the  inability  of  the  Kalriiau  liller  'o 
adequately  track  the  phase  On  as  shown  in  Fig.  10. 


Figure  11:  An  example  of  the  envelope-tracking  capability  of  the  Kalman  filter  when  the 
frequency  deviation  constant  d  was  varied  according  to  the  profile  in  Fig.  9. 

The  results  in  Figs.  10  and  11  suggest  that  an  optimum  value  of  d  exists  that  v-ill 
lead  to  accurate  tracking  of  the  interferer  and  minimum  mean-squared  interference  error; 
it  also  suggests  a  strategy  for  an  adaptive  algorithm.  An  architecture  that  illustrates  the 
concept  is  shown  in  Fig.  12,  which  is  based  on  finding  that  value  of  d  which  minimizes 
the  mean-squared  value  of  &t  the  input  to  the  PN  correlator.  This  concept  is  based 
on  the  use  of  the  interference  estimator  in  Fig.  7.  Figure  12  would  have  to  be  modified 
appropriately  by  incorporating  a  delay  of  D  samples  in  u„  and  the  PN  correlator  if  Fig.  S 
would  be  used;  the  results  presented  later  in  this  section  are  based  on  this  latter  approach. 


17 


Control 


PN 

CoireUtor 


Figure  12:  Adaptive  architecture,  which  is  based  on  minimizing  the  squared  error  of  the 
error  signal  Cni  while  varying  d. 

Minimizing  the  mean-squared  value  of  c„  is  equivalent  to  determining  the  mini¬ 
mum  of 


E{ei}  =  £{(s„-hn„  +  A.r} 

«  (38) 

where  is  the  residual  interference  which  has  been  defined  in  [5].  The  eissumptions 
which  were  made  so  that  the  terms  in  Eq.  (38)  could  be  de-coupled  are  also  contained  in 
[5]. 

Figure  13  conceptually  shows  the  steps  of  the  adaptive  process.  Figure  13(a) 
illustrates  an  example  of  the  d  profile  which  is  decreased  according  to  the  equation 

d,,  =  doexp"®*",  (39) 

where  do  is  the  initial  condition,  and  cq  controls  the  rate  of  change  of  d„.  This  part  of  the 
profile  forms  the  search  phase,  during  which  time  am  estimate  of  the  mean-squcired  value 
of  the  error  e„  at  the  input  to  the  PN  correlator  is  continually  being  calculated  using  the 
relation 
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Figure  13:  Dlustratioa  of  the  adaptive  Kalman  filter  process. 


where  is  a  number  close  to  1.0.  As  d  is  decreased  according  to  Eq.  (39),  a  scaled  version 
of  ^  (scaled  by  an  estimate  of  the  power  in  in  Fig.  12)  is  monitored.  This  scaled  lenn 
is  defined  as 

(41) 

where  is  an  estimate  of  the  power,  i.e., 

«l[  =  6«^+(l  ~6)«n  (42) 


where  is  a  number  close  to  1.0.  Equation  (42)  requires  a  reasonable  initial  condition  to 
counter  the  long  time  constants  resulting  from  those  values  of  (3  very  close  to  1.0.  One 
approach  that  is  used  in  the  results  to  follow  is  to  produce  a  ‘‘quick”  estimate  based  on 


the  relation 


(43) 


where  j  is  set  to  some  reasonable  value.  After  j  samples,  Eq.  (43)  is  replaced  by  Eq.  (42). 
Another  possibility  is  to  break  up  the  estimate  of  into,  say,  two  parts.  The  first  part 
could  correspond  to  one  value  of  ^3  that  would  produce  a  quick  estimate  of  and  the 
second  part  to  another  value  of  4a  closer  to  1.0,  using  the  result  from  the  first  part  as  an 
initial  condition. 

An  example  of  the  type  of  profile  created  by  Eq.  (41)  is  illustrated  in  Fig.  i3(b). 
The  search  for  the  minimum  commences  when  first  crosses  the  threshold  defined  by 
t/tresl  occurring  at  Ti  seconds.  The  value  of  d  (call  this  dgpt)  in  Fig.  13(a)  at  which  the 
minimum  of  occurs,  is  retained.  The  seuch  ceases  when  crosses  the  pre-defined 
threshold  denoted  by  thres2  at  Tj  seconds.  At  this  time  d  is  immediately  changed  to  the 
optimum  value  dgpt  determined  in  the  search  phase,  and  the  Kalman  filter  is  re-initialized. 

Two  examples  of  the  adaptive  algorithm  will  now  be  discussed  for  two  narrow- 
band  Gaussian  noise  interferers  of  bandwidths  Bi  =  0.01  Hz  and  0.05  Hz.  For  these 
examples,  Eh/No  =  12  dB,  the  interference-to-signal  ratio  is  20  dB,  and  L  =  20.  The  in¬ 
terference  estimator  was  the  one  in  Fig.  8.  The  low  pass  filter  had  7  taps  with  rectangular 
weighting  and  was  of  bandwidth  Bipf  —  0.20  Hz;  thus  the  delay  D  was  3  samples.  The 
phase  was  dso  smoothed  over  7  samples. 

A  few  words  must  be  mentioned  about  the  selection  of  specific  values  for  the 
constants  do  and  cq  in  Eq.  (39).  In  the  examples  to  be  considered,  do  has  been  set  to 
1.75  rad./sec./volt,  corresponding  to  a  reasonable  value  obtained  from  Fig.  14  ^  which 


^Thia  figure  shows  the  amount  of  interference  suppression  (i.e.,  S  =  10  log( Ai^/P,)  where  Pi  is  inter¬ 
ference  power)  for  various  values  of  d  and  interference  bandwidths,  Bt. 
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Figure  14:  Interference  suppression  level  5  as  a  function  d  for  several  interference  band- 
widtfis,  Bi,  ranging  from  0.01  Hz  to  0.05  Hz,  with  Bcpf  =  0.20  Hz. 
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was  extracted  from  [12]  for  the  interference  bandwidths  considered  therein.  The  value 
for  Co  is  selected  according  to  the  user’s  desired  maximum  change  in  d,  i.e.,  Ac/^ar. 
the  time-constant  determined  by  in  Eq.  (40)  which  estimates  the  power  in  the  error 
signal  e„;  for  the  examples  considered,  Ad„ax  was  chosen  to  be  0.10  rad. /sec. /volt,  and 
=  0.999  (a  time-constant  of  1000  seconds).  Therefore,  with  cq  =  0.000065  and  do  =  1.75 
r2ul./sec./volt,  d  was  guaranteed  not  to  change  by  more  than  0.10  rad./sec./volt  over  any 
1000  second  interval  during  the  search  phase.  This  strategy  provides  some  guidance  to  the 
user  as  to  the  rationale  in  selecting  values  for  these  constants,  bearing  in  mind,  however, 
that  several  trade-offs  involving  these  parameters  must  be  considered.  These  tradeoffs 
include  the  accuracy  in  estimating  ef  aind  uf,  speed  of  adaptation,  the  quality  of  the  r„„ 
curve  in  Fig.  13(b),  and  the  threshold  levels  thresl  and  thrt32. 

Figure  15  shows  the  profile  of  the  frequency  deviation  constant  that  resulted 
when  the  bandwidth  of  the  interferer  was  0.01  Hz.  The  steady-state  value  in  this  case 
was  dopt  —  0.30  rad./sec./volt  which  agrees  reasonably  well  with  the  results  in  Fig.  14. 
The  error  signal  Cn-o,  illustrated  in  Fig.  16,  shows  the  effect  of  the  changing  deviation 
constant  in  Fig.  15.  Ob.«erve  that,  as  d  becomes  too  small,  the  error  signal  starts  to 
increase,  indicating  less  interference  cancellation.  Tov/ards  the  end,  at  «  T  =  70, 000 
seconds,  when  fhres2  is  finally  crossed,  the  Kalman  filter  is  re-initialized  and  set  to  a 
value  of  dopty  resulting  in  a  lower  level  of  Cn-o.  The  mean-squared  value  of  the  reference 
signal  u«_o  is  illustrated  in  Fig.  17.  The  initial  condition  was  calculated  with  j  =  1000 
in  Eq.  (43).  After  1000  samples,  the  result  u*  was  inserted  into  Eq.  (42)  with  {2  =  0.9999 

The  mean-squared  value  of  en-o  was  calculated  using  Eq.  (40)  with  =  0.999 
and  scaled  by  uj  in  Fig.  17.  The  resulting  ratio  from  Eq.  (41)  is  shown  in  Fig.  18. 
The  two  thresholds  thresl  and  threa2,  were  set  to  0.05  (-13  dB)  and  0.07  (-11.6  dB) 
respectively;  these  values  were  determined  with  the  aid  of  Fig.  14  and  experimentally. 
Referring  to  Figs.  13  and  18,  thresl  is  crossed  at  Ti  »  4000  seconds  and  thres2  is  crossed 
at  T2  71000  seconds.  The  residual  interference  and  approximate  suppression  level  for 
this  case  aire  illustrated  in  Figs.  19  and  20.  The  latter  figure  shows  an  average  interference 
suppression  value  of  23  dB,  which  agrees  with  the  results  in  Fig.  14. 

Figures  21  to  25  illustrate  the  results  when  the  interference  bandwidth  has  been 
increased  to  Bi  =  0.05  Hz.  As  shown  in  Fig.  21,  the  final  “optimum”  value  for  d  is 
a;  0.82  rad./sec./volt.  This  larger  value  for  d  is  to  be  expected,  since  a  larger  bamdwidth 
is  required  for  the  Kalman  filter  to  track  the  faster  varying  phase.  Figure  22  illustrates 
the  power  estimated  in  the  delayed  input  signal  u„_o  which,  for  this  particular  example, 
provided  a  better  estimate  than  the  one  in  Fig.  17;  the  initial  condition  calculated  from 
Eq.  (43)  was  closer  to  the  desired  value.  The  ratio,  Vnitt  as  estimated  according  to 
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Eq.  (41),  is  illustrated  in  Fig.  23.  For  this  larger  bandwidth  inlerferer.  the  adaptive 
system  reacts  fcister  as  compared  to  the  case  where  B,  =  0.01  Hz.  the  of  which  was 
illustrated  in  Fig.  18.  The  residual  interference,  Ai„_o,  shown  in  Fig.  21.  compared  with 
the  results  in  Fig.  19,  reveals  a  significantly  higher  level  of  residual  interference.  Finally, 
Fig.  25  shows  the  average  interference  suppression  which  is  approximately  17  dB  after 
adaptation  is  completed. 

The  bit  error  performance  of  this  adaptive  algorithm  is  presented  next.  For  this 
example,  10,000  bits  were  transmitted  and  the  bit  errors  were  counted  over  time,  during 
and  after  the  adaptation  period.  The  conditions  were  the  same  as  those  used  in  the 
previous  examples  except  that  the  interferers  ranged  in  bandwidth  from  0.01  Hz  to  0.05 
Hz  in  steps  of  0.01  Hz.  The  bit  error  profiles  are  illustrated  in  Fig.  26.  The  profiles  consist 
of  four  sections  having  different  slopes,  which  correspond  to  the  high  and  low  levels  of  the 
residual  interference  results  in  Figs.  19  and  24.  Furthermore,  the  times  at  which  the 
adaptive  algorithm  switches  to  the  optimum  value  for  the  frequency  deviation  constant 
dopt  occur  in  the  pattern  as  expected:  as  Bi  increases  Tj  decreases. 

5.0  CONCLUDING  REMARKS 

This  report  has  presented  an  adaptive  approach  to  interference  suppression  using  a  Kalman 
filter.  The  Kalman  filter  is  also  referred  to  as  a  digital  phase-locked  loop,  since  the  sys¬ 
tem  structure  looks  like  an  analog  phase-locked  loop.  The  scheme  is  based  on  varying 
the  frequency  deviation  constant  d  of  the  Kaimaui  filter  over  time  while  holding  its  other 
parameters  fixed  and  finding  that  value  of  d  which  would  minimize  the  mean-squared 
error  at  the  excisor  output. 

The  results  show  that  this  approach  proves  to  be  a  valid  technique.  One  possible 
drawback  is  the  adaptation  time,  although  the  results  show  that  significant  interference 
rejection  is  achieved  for  much  of  the  adaptation  period.  Further  analysis  with  the  various 
system  parameters  would  ha.'e  to  be  conducted  in  order  to  establish  if  the  adaptation 
period  could  be  reduced,  bearing  in  mind,  of  course,  that  this  would  incur  penalties  in 
performance  because  of  the  requirement  to  balance  several  conflicting  tradeoffs. 

F’nally,  the  interference  which  was  considered  in  this  study  was  narrowband 
Gaussian  noise  whose  envelope  and  phase  both  vary  over  time.  For  constant  envelope 
interference  having  the  same  bandwidth  as  the  narrowband  Gaussian  noise,  one  would 
expect  better  performance  from  this  type  of  excisor. 
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Figure  15:  An  example  of  the  profile  of  the  frequency  deviation  constant  i 
2Mlaptive  algorithm,  Bi  =  0.01  Hz. 
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Figure  16:  An  example  of  the  error  signal  e„_D  for  the  d  profile  in  Fig.  15  used  in  the 
2ulaptive  algorithm,  with  Bi  =  0.01  Hz. 
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Figure  17:  An  example  of  the  estimated  power  in  Un-o  C2jculated  from  Eq.  (42)  for  the 
d  profile  in  Fig.  15  used  in  the  adaptive  algorithm,  with  Bi  =  0.01  Hz. 
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Figure  19:  An  example  of  the  residual  interference  Ai„_D  =  in-D  —  in-D  before  the 
despreader  in  the  2tdaptive  2dgorithm,  with  Bi  =  0.01  Hz. 
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Figure  20:  An  example  illustrating  the  degree  of  interference  suppression 

101og(AiJ_|j/P,),  where  P,  =  100  is  the  power  in  the  interferer  of  bandwidth  Bi  =  0.01 


Figure  21:  An  example  of  the  profile  of  the  frequency  deviation  constant  d  used  in  the 
adaptive  2dgorithm,  Bi  =  0.05  Hz. 
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Figure  22:  An  example  of  the  estimated  power  in  Un-D  calculated  from  Eq.  (42)  for  the 
d  profile  in  Fig.  21  used  in  the  adaptive  algorithm,  with  B,  =  0.05  Hz. 
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Figure  23:  An  example  of  the  ratio  of  rn,e  =  ^n-ol^n-o  profile  in  Fig.  21  used 

in  the  adaptive  algorithm,  with  B,  =  0.05  Hz. 
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Figure  25:  An  example  illustrating  the  degree  of  interference  suppression 

101og(AtJ~^/Pi),  where  Pi  =  100  is  the  power  in  the  interferer  of  bandwidth  B,  =  0.05 


Hz. 
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Figure  26:  Bit  error  performance  of  the  adaptive  scheme  for  several  interference  band- 
widths  and  EifNo  =  12  dB  using  the  interference  estimator  in  Fig.  8. 
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